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The  boundary  conditions  (BCs)  implemented  in  the  Pacific  West  Coast  (PWC)  version  of  the  Princeton  Ocean  Model  (POM)  as 
developed  at  the  Naval  Research  Laboratory  (NRL)  are  described.  BCs  for  the  barotropic  velocities  along  the  open  boundaries  are  a 
Flather  form  for  the  normal  component  and  an  advective  form  for  the  tangential  component  with  matching  to  externally  supplied 
transports.  Baroclinic  velocities  use  a  radiational  BC  for  the  normal  component,  and  an  advectional  BC  for  the  tangential  component. 
Sea  surface  height,  temperature,  and  salinity  use  radiational  or  advective  BCs.  Simple  relaxation  techniques  are  used  for  the  assimi¬ 
lation  of  sea  surface  temperature  (SST)  and  sea  surface  salinity  (SSS).  The  interior  state  of  the  ocean  is  maintained  by  relaxation  to 
a  monthly  climatology.  A  zero  gradient  boundary  condition  is  applied  for  the  vertical  velocity,  a  clamped  condition  is  applied  to  the 
turbulent  kinetic  energy  and  turbulent  length  scale.  The  explicit  treatment  of  these  BCs  is  given  here  to  serve  as  a  reference  for  PWC 
POM  users. 
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Boundary  Conditions  in  the  Pacific  West  Coast 
Princeton  Ocean  Model  of  CoBALT 


1.  INTRODUCTION 

An  ocean  general  circulation  model  (OGCM)  for  the  Pacific  West  Coast  (PWC)  was  developed 
by  NRL  in  an  effort  to  meet  Navy  operational  needs  for  this  region  [Clancy  et  al.,  1996].  The 
OGCM  was  based  on  the  POM  with  extensive  modifications  to  include  one-way  coupling  from  the 
NRL  Layered  Ocean  Model  (NLOM)  [ Ko ,  1996a;  Ko  et  al.,  1996b;  Ko  et  al.,  1996c;  Metzger  et 
al,  1998].  The  PWC  POM  as  originally  developed  by  Ko  demonstrated  sufficient  predictive  skill 
that  it  was  delivered  to  Fleet  Numerical  Meteorology  and  Oceanography  Center  for  evaluation  as 
an  operational  product  for  use  by  the  fleet.  Since  then  this  version  has  been  coupled  to  a  regional 
weather  forecast  model  [Allard  et  al.,  1996],  been  turned  into  a  nowcast /forecast  model  with  data 
assimilation  [Ko  et  al.,  1999a],  and  coupled  to  an  OGCM  of  the  North  Pacific  Ocean  [Ko  et  al., 
1999b].  A  close  variant  of  the  PWC  POM  is  presently  being  used  for  studies  as  part  of  the  Coupled 
Biophysical-Dynamics  Across  the  Littoral  Transition  (CoBALT)  project  at  NRL  [Haidvogel  et  al., 
2000;  Kindle  et  al.,  1999;  Righi  et  al.,  1999;  Shulman  et  al.,  1999a]  The  CoBALT  version  of  the 
PWC  POM  differs  from  the  original  due  to  modifications  made  to  the  BCs  for  the  external  mode  as 
devised  by  Shulman  et  al.  [1999b]  In  late  1999  the  real-time  simulation  at  FNMOC  was  modified 
to  be  consistent  with  the  CoBALT  version  of  the  PWC  POM.  Given  the  wide  applicability  of 
this  OGCM,  and  the  uniqueness  of  the  boundary  conditions  that  embed  the  coastal  model  within 
the  global  NLOM,  we  document  in  this  report  the  details  of  the  implementation  of  the  boundary 
conditions  for  the  benefit  of  the  user  community.  The  BCs  described  are  those  specifically  used 
in  the  PWC  POM  of  the  CoBALT  project,  which  we  will  hereinafter  referred  to  as  just  the  PWC 
POM. 

The  open  boundaries  along  the  south,  north,  and  western  edges  of  the  PWC  POM  domain  are 
forced  by  fields  provided  from  the  global  NRL  Layered  Ocean  Model  (NLOM).  The  PWC  POM  is 
bounded  on  the  eastern  edge  of  the  model  domain  by  the  Pacific  west  coast  of  the  United  States  and 
includes  temperature  and  salinity  forcing  as  dictated  by  river  inflows.  The  model  also  assimilates 
sea  surface  temperature  (SST)  from  daily  predictions  of  Multi-Channel  Sea  Surface  Temperature 
(MCSST)  and  sea  surface  salinity  (SSS)  from  the  monthly  climatology  of  the  World  Ocean  Atlas 
1994  [Levitus  et  al.,  1994;  Levitus  and  Boyer,  1994]  (hereinafter  referred  to  as  Levitus  data).  The 
state  of  the  interior  ocean  is  maintained  via  a  relaxation  to  monthly  Levitus  temperature  and 
salinity. 

Section  2  presents  the  BCs  as  implemented  for  elevation.  Sections  3  and  4  describe  the  BCs 
for  the  barotropic  and  baroclinic  velocities,  respectively.  Section  5  explains  the  temperature  and 
salinity  BCs,  the  assimilation  of  MCSST  and  Levitus  SSS,  and  the  relaxation  of  the  interior  tem¬ 
perature  and  salinity  to  climatology.  Section  6  presents  the  BCs  for  the  vertical  velocity.  Section  7 
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gives  the  BCs  for  the  turbulent  kinetic  energy  and  the  turbulence  length  scale.  Finally,  in  section  8 
we  provide  a  short  summary  of  the  BCs. 

2.  ELEVATION 


For  the  elevation  r?  a  Sommerfeld  radiation  BC  is  applied 


dn  dy 
dt  +Cbtdn 


=  0, 


(1) 


where  n  is  the  unit  outward  normal  to  the  boundary,  and  Cbt  is  the  barotropic  phase  speed.  For 
the  PWC  region  n  =  —y  for  the  southern  boundary,  n  —  y  for  the  northern  boundary,  and  n  =  —  x 
for  the  western  boundary.  The  phase  speed  is  given  by 

Cbt  =  {gH)ll\  (2) 


where  g  is  the  gravitational  acceleration  and  H  is  the  total  depth  of  the  ocean. 

In  the  FORTRAN  implementation  of  this  BC  (see  BCOND  subroutine  in  POM)  a  1:2:1  filter 
is  applied  to  the  elevations  at  the  boundary  before  the  BC  is  applied  in  order  to  reduce  noise 

Vx(iJ)  =  \  (vii  ~  1  ,j)  +  2 y(i,j)+r](i  +  l,j))  (3) 

Vy(i,j)  =  \  (v(h3  ~  1)  +2 r](i,j)  +r](i,j  +  1)) ,  (4) 

where  fjx  and  fjy  denote  filtering  in  the  zonal  and  meridional  directions,  respectively.  The  elevation 
at  the  boundaries  is  updated  for  the  forward  time  step  t  +  At  via 

fj(t  + At)  =  y(t)  -  Cbt^^-AtE  (5) 

where  Arj  represents  the  finite  difference  of  the  elevation  and  A tE  is  the  external  time  step.  This 
boundary  condition  is  applied  to  the  southern,  northern,  and  western  boundaries. 


As  an  example,  the  FORTRAN  code  for  the  elevation  BC  at  the  southern  boundary  is  given 
below. 


DO  I  =  2,  NBS 

CEES  (I)  =  SQRT(  (H(I,1)+H(I,2))*.5*GRAV  ) 

&  *  DTE  /  (  (DY (I , 1)+DY (I ,2))* .5  ) 

ELF(I.l)  =  CEES (I)  * ( .25*EL(I-1 ,2)+.5*EL(I ,2)+. 25*EL(I+1 ,2) ) 

&  +(1-CEES(I))*(.25*EL(I-1,1)+.5*EL(I,1)+.25*EL(I+1,1)) 

END  DO 

The  FORTRAN  symbols  are  defined  in  Appendix  A.  Note  that  H  and  DY  are  averaged  to  colocate 
them  with  V  on  the  C-grid  (Appendix  B). 

! 

The  elevation  at  the  SW  and  NW  corners  of  the  model  grid  are  given  special  treatment,  as 
they  are  averaged  with  their  boundary  neighbours. 
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ELF(1 , 1  )  =  (ELF(1 ,2)+EL(2, l))/2. 

ELF(l,jm)  =  (ELF(l,JMMl)+EL(2,JM))/2. 

This  is  done  for  numerical  safety  in  the  event  these  corner  locations  are  accessed  elsewhere  in  the 
model. 

For  numerical  safety  the  elevations  are  also  specified  at  the  river  locations  along  the  coast. 
Their  values  are  assigned  to  those  of  the  adjacent  interior  grid  point  where  they  are  computed  by 
the  POM  because  there  is  no  suitable  information  on  river  elevations  at  ocean  boundaries.  This  is 
coded  into  the  POM  as 

DO  I  =  1,  NRV 

ELF(IRV(I) ,  JRV(I))  =  ELF(IRV(I)-1,JRV(I)) 

END  DO 

where  IRV  and  JRV  are  grid  locations  of  the  rivers. 

3.  BAROTROPIC  VELOCITIES 


For  the  barotropic  velocity  V  =  (U,  V)  different  BCs  are  applied  for  the  normal  and  tangential 
components.  For  the  normal  component  Vn  a  radiation  BC  that  was  used  for  barotropic  two- 
dimensional  tide  and  storm  surge  modeling  [Flather,  1976]  is  applied 

Vrn  =  ^„°+(|:)1/2(r?-7?o),  (6) 


where  V®  is  the  NLOM  normal  velocity  at  the  open  boundary,  77  is  the  model  sea  surface  elevation 
calculated  from  the  continuity  equation,  7/0  is  the  NLOM  elevation  at  the  open  boundary.  Note  that 
Flather  [1976]  considered  both  meteorological  (Vmet  and  r;rnet)  and  tidal  forcing  (Vtide  and  7?tide)  in 
the  BC 


Vn  —  ^met  "I"  I'tide  — ' 


(V  rlmet  —  f/tide)- 


(7) 


For  the  tangential  component  Vt  an  advectional  BC  is  used  to  allow  one  way  external  forcing 
of  the  PWC  POM.  The  vertically  averaged  tangential  component  of  velocity  from  the  NLOM  is 
advected  into  the  model  domain  in  the  case  of  inflow,  and  the  internal  vertically  averaged  tangential 
component  of  velocity  is  advected  to  the  open  boundary  in  the  case  of  outflow.  This  leads  to  a  BC 
of  the  form 


dvt 

dt 


(8) 


where  here 


(Vt-Vg)/ Axn  vn>0 

On  \  (Vi  -  Vt)/Axn  Vn  <  0  W 

{ 

f 

with  Vg  =  Vto/HAxn  the  tangential  barotropic  velocity  at  the  open  boundary  obtained  from 
the  NLOM  transport  Vjo,  V{  the  tangential  barotropic  velocity  at  one  grid  point  inside  the  open 
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boundary,  and  Axn  the  grid  spacing  in  the  direction  normal  to  the  boundary.  The  Flather  BC  is 
applied  to  V  at  the  southern  and  northern  boundaries  and  to  U  at  the  western  boundary.  The 
advectional  BC  is  applied  to  U  at  the  southern  and  northern  boundaries  and  to  V  at  the  western 
boundary. 

As  an  example,  the  FORTRAN  code  for  the  normal  barotropic  BC  at  the  southern  boundary 
is 

DO  I  =  2,  NBS 

DEPTH=(  H(I,1)+  H(I ,2)  )/2. 

DV  =(  DX(I,1)+DX(I,2)  )/2 . 

C0VRHS(I)=  SQRT(  GRAV/H(I,2)  )*DEPTH*DV 

VAF(I,2)  =  (VAGS(I)*RAMP  - 

&  COVRHS(I) * (  EL ( I , 2 ) -ELGS ( I ) *RAMP  )  )/(DEPTH*DV) 

ENDDO 

while  that  for  the  tangential  barotropic  BC  is 
DO  I  =  2,  NBS 

VMID  =  (  VA(I-1 ,2)+VA(I ,2)  )/2. 

DEPTH=  (  H(I-1 , 1)+H(I , 1)  )/2 . 

DV  =  (  DY(I-1 , 1) +DY(I , 1)  )/2. 

UAG  =  RAMP*UAGS ( I ) / (DEPTH*DV ) 

AEUS(I)  =  DTE/ (DY(I-1 , 1)+DY(I , 1)+DY(I-1 ,2)+DY(I , 2) ) *2 . 

UAF(I , 1)=  UA(I,1)  -  AEUS(I)  * 

&  (  (VMID+ABS (VMID) )  *  (UA ( I , 1 ) -UAG  ) 

&  + (VMID-ABS (VMID) )  *  (UA(I ,2)-UA(I , 1) )  ) 

END  DO 

Note  that  H  in  DEPTH,  V  in  VMID,  and  Axn  in  DV  are  averaged  over  2  grid  points  to  colocate 
at  the  corresponding  U  and  V  points  on  the  C-grid  (Appendix  B)  the  phase  speed  (g/H)1/2  and 
normal  velocity  Vn ,  respectively.  Similarly,  the  grid  space  averaging  of  Ay  in  AEUS  is  to  properly 
account  for  the  tangential  velocity  derivative  dVt/dn. 

After  the  BCs  have  been  applied  to  the  barotropic  velocities  the  values  of  V  at  the  first  row 
of  the  southern  boundary  and  the  values  of  U  at  the  first  column  of  the  western  boundary  are 
assigned  the  values  of  their  adjacent  grid  points  for  numerical  safety. 

DO  I  =  1,  NBS 

VAF(I , 1)  =  VAF(I,2) 

END  DO 

DO  J  =  1,  NBW 

UAF(1 , J)  =  UAF(2 , J) 

END  DO 

River  runoff  is  included  as  a  BC  for  the  barotropic  velocities.  Only  U  is  specified  because  all 
the  rivers  flow  from  east  to  west.  It  is  coded  as 


DO  IJ  =  1,  NRV 
I  =  IRV(IJ) 

J  =  JRV(IJ) 
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DEPTH  =  (H(I-1, J)+ELF(I-1, J)+H(I, J)+ELF(I, J))/2.0 
DV  =  (DY(I-1, J)+DY(I, J))/2.0 
UAF(I.J)  =  RAMP*UARV(I J) / (DEPTH*DV) 

END  DO 

where  UARV  is  the  river  transport. 

4.  BAROCLINIC  VELOCITIES 


The  baroclinic  velocity  v  —  (u,v)  has  a  radiational  BC  applied  for  the  normal  component  vn 
that  is  identical  in  form  to  that  used  for  the  elevation 


dvn  ,r  dvn 

IT  +  c"‘~, *T  =  °' 


(10) 


but  where  the  phase  speed  is  assigned  as  a  fraction  of  the  barotropic  phase  speed 

Cbc  =  0.005  {gH)1/2.  (n) 

This  explicit  relation  was  determined  through  numerical  tests  to  yield  adequate  results.  For  the 
tangential  component  vt ,  an  advectional  BC  is  applied  that  is  identical  in  form  with  that  used  for 
the  barotropic  tangential  velocity 


dvt  dvt 


The  normal  gradient  of  the  velocity  is  given  here  by 

_  f  (vt  ~  vg)/Axn  vn>0 
dn  \  (vi  -  vt)/Axn  vn  <  0 


(12) 


(13) 


where  Vi  is  the  tangential  barotropic  velocity  at  one  grid  point  inside  the  open  boundary,  and  vg  =  0 
because  there  are  no  vertical  profiles  of  tangential  velocity  available  at  the  open  boundary.  The 
radiational  BC  is  applied  to  v  at  the  southern  and  northern  boundaries  and  to  u  at  the  western 
boundary.  The  advectional  BC  is  applied  to  u  at  the  southern  and  northern  boundaries  and  to  v 
at  the  western  boundary. 


The  FORTRAN  implementation  of  the  BC  for  the  normal  baroclinic  velocity  at  the  southern 
and  northern  boundaries  are  given  below  as  examples. 

ck  **  South 

DO  I  =  2,  NBS 

CIVS(I)  =  SQRT (H (I , 2) *GRAV*5 . E-3) *DTI/DY(I , 2) 

VF(I,2,K)  =  CIVS(I)  *(.25*V(I-1,3,K)+.5*V(I,3,K)+.25*V(I+1,3,K)) 

&  +(1“CIVS(I))*( .25*V (1-1,2, K)+.5*V( 1,2, K) +.25* V( 1+1,2, K)) 

END  DO 

ck  **  North 

DO  I  =  2,  NBN  , 

CIVN  (I)  =  SQRT (H ( I , JMM1 ) *GRAV*5 . E-3 ) *DTI/DY ( I , JMM1 ) 

VF(I , JM,K)  =  CIVN (I)  *  (.25*V(I-1, JMM1,K)+.5*V(I, JMMl.K) 
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&  +.25*V(I+1 , JMM1 ,K) ) 

&  +(1-CIVN(I))*(.25*V(I-1, JM,K)+.5*V(I, JM,K) 

&  +.25*V(I+1 , JM,K) ) 

END  DO 


The  corresponding  coding  for  the  tangential  baroclinic  velocity  at  the  southern  and  northern  bound¬ 
aries  are  as  below. 


ck  **  South 

DO  I  =  2,  NBS 

AIUS(I)  =  2 . *DTI/ (DY (I— 1 , 1)+DY(I , 1)+DY(I-1 ,2)+DY(I ,2) ) 
VMID  =(  V(I,2,K)+V(I-1,2,K)  )/2. 

UF(I , 1 ,K)=U(I , 1 ,K)  -  AIUS(I)  * 

&  (  (VMID+ABS(VMID))*(U(I,1,K)-0.E0  ) 

&  +(VMID-ABS(VMID))*(U(I,2,K)-U(I,1,K))  ) 

END  DO 


ck  **  North 

DO  I  =  2,  NBN 

VMID=.5E0*(V(I,JM,K)+V(I-1,JM,K)) 

UF(I , JM,K)=U(I , JM,K)  -  AIUN(I)  * 

&  (  (VMID+ABS(VMID) ) * (U(I , JM ,K) -U(I , JMM1 ,K) ) 

&  +(VMID-ABS (VMID) ) * (0 . EO  -U(I,JM,K)  )  ) 

END  DO 

Note  the  2  interior  points  closest  to  the  domain  boundary  are  used  for  the  normal  BCs  at  the 
southern  and  western  boundaries  because  of  the  way  the  C  grid  is  defined  for  the  POM  (see 
Appendix  B). 


As  in  the  case  of  the  barotropic  velocities  the  values  of  v  at  the  first  row  of  the  southern 
boundary  and  values  of  u  at  the  first  column  of  the  western  boundary  are  assigned  the  values  of 
their  adjacent  grid  points  for  numerical  safety. 

DO  I  =  1,  IM 

VF(I,1,K)  =  VF(I,2,K) 

END  DO 

DO  J  =  1,  JM 

UF(1,J,K)  =  UF(2,J,K) 

END  DO 

There  is  no  suitable  information  on  vertical  profiles  of  river  inflows  to  use  for  the  baroclinic 
velocities.  The  values  of  u  at  the  river  locations  are  therefore  assigned  to  those  of  the  adjacent 
interior  grid  point  where  they  are  computed  by  the  POM.  This  is  coded  as  below. 

DO  I  =  1,  NRV  i 

UF(IRV(I),JRV(I) ,K)  =  UF(IRV(I)-1,JRV(I),K) 

END  DO 
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5.  TEMPERATURE  AND  SALINITY 


5.1  Water  Column 


An  advectional  BC  is  applied  along  the  open  boundary  for  the  temperature  T  and  salinity 
S  throughout  the  water  column.  It  is  identical  in  form  to  that  used  in  (12)  for  the  tangential 
baroclinic  velocity 


dT  &r 

dt  Vn  dn 


0 


dS  dS  n 

m+Vnd^  =  0 


(14) 

(15) 


The  normal  gradients  of  the  temperature  and  salinity  are  given  here  by 

dT  =  f  (T  -  T0)/Axn  vn>0 
dn  {  ( Ti  -  T)/Axn  vn  <  0 

dS_  f  (S-S0)/Axn  vn>0 
dn  (  (Si  —  S)/Axn  vn  <  0 


(16) 

(17) 


where  To  and  So  here  are  the  climatological  temperature  and  salinity  at  the  boundary,  and  T)  and 
Si  the  temperature  and  salinity  at  one  grid  point  inside  the  open  boundary,  respectively.  The 
advectional  BC  is  applied  to  T  and  S  at  the  southern,  northern  and  western  boundaries. 


The  FORTRAN  implementation  of  these  BCs  at  the  southern  boundary  is  as  below. 

DO  I  =  2,  NBS 

ATSS(i)  =  DTI/ (DY(I , 1)+DY (1,2)) 

UF(I , 1 ,K)  =  T(I , 1 ,K)  -  ATSS(I)  * 

&  (  (  V (1 ,2 ,K)+ABS(V (1 ,2  ,K)  )  )*(T(1 , 1 ,K)-TBS (I ,K) ) 

&  +(  V ( I , 2 , K) -ABS (V ( I , 2 , K)  )  )*(T(I,2,K)-T(I,1,K))  ) 

VF(I,1,K)  =  S(I , 1 ,K)  -  ATSS(I)  * 

&  (  (  V (I ,2 ,K)+ABS(V(I ,2 ,K)  )  )*(S(I,1,K)-SBS(I,K)) 

&  +(  V(I,2,K)-ABS(V(I,2,K)  )  )*(S(I,2,K)-S(I,1,K))  ) 

END  DO 

Note  that  the  temporary  variables  UF  and  VF  are  used  to  contain  T  and  S  on  entry  into  the 
BCOND  subroutine. 


As  in  the  case  of  the  elevation,  the  SW  and  NW  corners  of  the  model  grid  are  averaged  with 
their  boundary  neighbours. 

UF(1 ,  1 ,K)  =  (UF(2,l,K)+UF(l,2,K))/2. 

VF(1 ,  1 ,K)  =  (VF(2,l,K)+VF(l,2,K))/2. 

UF(1,JM,K)  =  (UF(2,JM,K)+UF(l,JMMl,K))/2. 

VF(1 , JM,K)  =  (VF(2, JM,K)+VF(1, JMMl,K))/2. 

The  temperature  and  salinity  at  the  mouths  of  rivers  is  specified  in  the  PWC  POM.  For 
temperature  the  value  is  specified  to  be  that  given  by  river  data,  with  h  blending  of  the  ocean  and 
river  temperatures  being  used  during  the  ramp  up  phase  of  the  model.  For  salinity  the  river  is 
considered  to  be  pure  fresh  water.  This  BC  is  coded  as 
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DO  I  =  1,  NRV 
DO  K  =  1,  KBM1 

UF(IRV(I) , JRV(I) ,K) 


VF(IRV(I) ,JRV(I) ,K) 
END  DO 
END  DO 


UF(IRV(I) ,JRV(I) ,K) 

+  (TRV(K.I)  -  UF(IRV(I) , JRV(I) ,K)) 
*  RAMP 

VF(IRV(I)  ,JRVd)  ,K)*(1.-RAMP) 


where  TRY  is  the  temperature  of  the  river. 


5.2  Sea  Surface  Temperature  and  Salinity 


The  surface  temperature  and  salinity  used  in  the  POM  are  the  values  of  the  profiles  in  the 
previous  subsection  with  a  relaxation  to  externally  supplied  fields  of  SST  and  SSS.  In  the  original 
implementation  a  spatially  filtered  SST  field  ( SSTj )  is  first  produced  from  the  surface  values  of  the 
profiles  using  a  gaussian  weight  distribution  having  an  attenuation  length  of  5  grid  points.  This 
was  done  because  the  original  MCSST  used  was  created  from  5  day  composites  of  the  satellite 
images  and  was  rather  noisy.  The  PWC  POM  of  CoBALT  uses  daily  MCSST  from  the  Modular 
Ocean  Data  Assimilation  System  which  do  not  have  such  problems.  No  spatial  filtering  is  therefore 
applied  in  the  PWC  POM  when  using  daily  MCSST  ( SSTj  =  MCSST).  No  spatial  filtering  is 
applied  to  the  SSS  values  of  the  profiles.  The  sea  surface  temperature  and  salinity  fields  (SSTq  and 
SSSq)  are  then  assimilated  into  the  model  via  a  relaxation  method 


dSST  .  „  , 

dt  =  ^ t(SST0  —  SSTj) 

(18) 

8SSS 

=  As(SSS0  -  SSS). 

(19) 

In  the  original  implementation  of  the  PWC  POM  time  scales  of  AJ1  =  5d  and  A^1  =  30 d  were 
used.  When  the  PWC  POM  is  forced  with  daily  MCSST  the  value  A^1  =  Id  is  used  instead. 
(Note  that  a  two  timesteps  leapfrog  scheme  is  used  in  the  finite  difference  forms  of  the  POM.)  The 
SSTq  and  SSSq  fields  are  linearly  interpolated  to  the  model  time  step.  Note  that  the  relaxation  is 
applied  independent  of  the  model  time  step.  The  PWC  POM  presently  assimilates  daily  MCSST 
and  monthly  Levitus  SSS. 


5.3  Interior  Relaxation  to  Climatology 


To  prevent  model  drift  away  from  the  seasonal  cycle  the  subsurface  temperature  and  salinity 
are  relaxed  towards  the  Levitus  monthly  mean  temperature  (T)  and  salinity  ( S ) 


f  =  ast(t-D 

(20) 

—  =  A st(S  -  S). 

(21) 

The  time  scale  for  the  relaxation  is  given  by 


(22) 
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where  tsr  =  250  d,  zq  =  500m,  and  z  <  0.  The  250  day  time  scale  was  found  to  be  adequate  for 
having  the  interior  temperature  and  salinity  follow  the  seasonal  cycle.  The  time  scale  is  exponen¬ 
tially  weighted  with  depth  to  ensure  T  and  S  at  the  surface  only  relaxes  to  the  externally  supplied 
SST  and  SSS  at  the  surface.  A  check  is  subsequently  made  on  the  resulting  values  of  salinity  and 
values  of  S  <  0  are  set  to  zero. 

6.  VERTICAL  VELOCITY 

A  zero  gradient  BC  in  the  horizontal  is  applied  to  the  vertical  velocity  w  along  all  3  open 
boundaries  of  the  model  domain 


Vhw  =  0  (23) 

where  V#  is  the  horizontal  gradient.  Vertical  velocities  at  river  locations  are  assigned  to  those  of 
the  adjacent  interior  grid  point  where  they  are  computed  by  the  POM. 

These  BCs  are  coded  within  the  model  as  below. 

DO  245  K=1 ,KBM1 
DO  246  J=1 , JM 

246  W(l, J,K)=W(2, J,K) 

DO  247  1=2, IM 

W(I,JM,K)=W(I, JMM1 ,K) 

247  W(I,1,K)=W(I,2,K) 

245  CONTINUE 

C 

C  Rivers 

DO  I  =  1,  NRV 
DO  K  =  1,  KBM1 

W(IRV(I) , JRV(I) ,K)  =  W(IRV(I)-1, JRV(I) ,K) 

END  DO 
END  DO 

7.  TURBULENT  KINETIC  ENERGY  AND  LENGTH  SCALE 

A  clamped  BC  is  applied  to  twice  the  turbulent  kinetic  energy  q 2  and  its  product  with  the 
turbulence  length  scale  l  along  the  open  boundaries.  Values  of  q2  and  q2l  at  rivers  are  assigned 
those  of  the  adjacent  interior  grid  point. 

The  FORTRAN  coding  of  these  BCs  is  given  below. 

DO  300  K=1 ,KB 
DO  295  J=l, JM 

UF(IM, J,K)=1 .E-10 
VF(IM, J,K)=1.E-10 
DO  296  1=1, IM 

UF(I , JM,K)=1 .E-10  { 

VF(I , JM,K)=1 .E-10 
UF(I,1,K)=1.E-10 


295 
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296  VF(I,1,K)=1.E-10 

C 

C  Rivers 

DO  I  =  1,  NRV 

UF(IRV(I) , JRV(I) ,K)  =  UF(IRV(I)-1,JRV(I),K) 

VF(IRV(I),JRV(I),K)  =  VF(IRV(I)-1,JRV(I),K) 

END  DO 
300  CONTINUE 

Note  that  the  temporary  variables  UF  and  VF  are  used  to  hold  q 2  and  q2l  on  entry  into  the  BCOND 
subroutine. 

8.  SUMMARY 

We  have  presented  here  the  explicit  details  of  how  externally  supplied  forcing  has  been  im¬ 
plemented  for  one  way  coupling  along  the  open  and  coastal  boundaries  of  the  PWC  POM.  These 
forcing  fields  derive  from  a  variety  of  sources  that  include  a  coarser  resolution  OGCM,  satellite 
derived  fields,  and  river  inflows.  A  variety  of  formulations  are  implemented  in  the  PWC  POM  that 
include  advectional,  clamped,  Flather,  and  radiational  BCs,  as  is  appropriate  for  each  variable. 
Most  of  these  BCs  have  been  implemented  in  some  form  in  the  Navy  Coordinate  Ocean  Model 
(NCOM)  and  it  is  hoped  that  this  report  will  be  of  value  to  NCOM  users  as  well. 
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Appendix  A 


FORTRAN  SYMBOLS 


The  FORTRAN  variable  names  used  in  the  POM  that  are  cited  within  this  report  are  listed 
below  for  ease  of  reference.  The  variables  are  listed  as  given  in  the  POM  Users  Guide  [Mellor, 
1998]  and  the  reader  is  referred  to  the  latter  for  a  more  complete  list.  The  variables  are  given 
followed  by  their  corresponding  analytical  symbols  in  parentheses  and  a  brief  explanation.  Not 
explicitly  tabulated  are  the  suffixes  B,  blank  and  F  which  are  appended  to  many  of  the  vari¬ 
ables  to  denote  the  previous  (n  —  1),  latest  (n),  and  forward  (n  +  1)  time  levels,  respectively. 


Indices 

U  (id) 

horizontal  grid  indices 

IM,JM 

outer  limits  of  I  and  J 

IMM1,JMM1 

IM-1  and  JM-1 

IRV,JRV 

grid  indices  of  river  locations 

K  (k) 

vertical  grid  index;  K=1  at  the  top  and  K=KB  at  the 
bottom 

NBS,  NBN,  NBW 

the  total  number  of  ocean  grid  points  along  the  S,  N,  and  W 
boundaries 

Constants 

DTE  (A tE) 

external  mode  time  step  (s) 

DTI  (At/) 

internal  mode  time  step  (s) 

GRAV  (y) 

acceleration  due  to  gravity  (ms-2) 

RAMP  (A) 

amplitude  during  spin  up  from  rest  (0  <  A  <  1) 

One-dimensional  Arrays 

ELGS  (ijo) 

externally  elevation  at  open  boundary  (m) 

UAGS,  VAGS  (U0,V0) 

externally  supplied  transport  at  open  boundary  (ms-1) 

Two-dimensional  Arrays 

DX  (As) 

longitude  grid  spacing  (m) 

DY  (Ay) 

latitude  grid  spacing  (m) 

EL  (V) 

surface  elevation  as  used  in  the  external  mode  (m) 

H  (H) 

the  bottom  depth  (m) 

TRV 

the  temperature  of  the  river  (°C  ) 

UA,  VA  ( U ,  V) 

vertical  mean  of  v  (ms-1) 

Three-dimensional  Arrays 

u,  V(ti,t>) 

horizontal  velocities  (ms  x) 

W  (u) 

sigma  coordinate  vertical  velocity  (ms-1) 
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Appendix  B 


PRINCETON  OCEAN  MODEL  C-GRID 


The  POM  uses  the  C-grid  shown  in  Figure  Bl.  The  locations  of  the  velocity  components  are 
indicated  by  U  and  V,  and  the  mass  point  locations  by  the  elevation  E.  The  temperature  and 
salinity  are  at  the  mass  points  and  are  therefore  coincident  with  E.  The  diagram  shown  is  for  the 
case  where  velocities  are  the  primary  boundary  conditions  together  with  temperature,  salinity,  and 
elevation.  The  elevation  itself  is  rather  unimportant  and  is  substituted  from  an  adjacent  interior 
point. 

The  symbols  appearing  within  the  dotted  box  correspond  to  the  interior  (non-boundary)  grid 
points.  In  general  only  those  variables  in  the  interior  are  computed  and  variables  at  the  open 
boundary  have  to  be  specified.  All  interpolations  are  centered  in  space  except  those  at  lateral  open 
boundaries  where  an  upstream  scheme  is  usually  used.  The  BC  symbols  indicate  the  boundaries 
for  the  various  variables  on  the  C-grid.  The  rows/columns  indicated  by  NU  are  not  used  in  the 
model  calculations  as  are  variables  marked  by  an  asterix  (*).  Note  that  adjacent  interior  values 
may  be  filled  in  at  these  latter  points  for  more  attractive  output  of  the  fields. 
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Fig.  B1  -  The  C-grid  Used  by  the  PVVC  POM. 
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